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Abstract 

The multiplicity fluctuations are studied in the van der Waals excluded volume hadron-resonance 
gas model. The calculations are done in the grand canonical ensemble within the Boltzmann statistics 
approximation. The scaled variances for positive, negative and all charged hadrons are calculated along 
the chemical freeze-out line of nucleus- nucleus collisions at different collision energies. The multiplicity 
fluctuations are found to be suppressed in the van der Waals gas. The numerical calculations are 
presented for two values of hard-core hadron radius, r = 0.3 fm and 0.5 fm, as well as for the upper 
limit of the excluded volume suppression effects. 
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I. INTRODUCTION 



During last decades statistical models of strong interactions have served as an important tool 
to investigate high energy nuclear collisions. The main obiect of the study have been the 

n n n 

mean multiplicities of produced hadrons (see e.g. Refs. 12|, |3|). Only recently, due to a 
rapid development of experimental techniques, first measurements of fluctuations of particle 
multiplicities and transverse momenta {sl were performed. The study of event-by-event 
fluctuations in nucleus-nucleus collisions (see e.g., reviews ^) is motivated by expectations of 
anomalies in the vicinity of the onset of deconfinement |3] and in the case when the expanding 
system goes through the transition line between the quark-gluon plasma and the hadron gas 18] . 
Furthermore, it is expected that a critical point of strongly interacting matter may be signaled 
by a characteristic pattern in fluctuations [9^. 

A convenient measure of the particle number fluctuations is the scaled variance, u = 
Var{N)/{N), defined in terms of the mean value, (N), and variance, Var{N) = {N"^) — (N)"^. 
The scaled variance equals to 1 for the ideal Boltzmann gas in the grand canonical ensemble 
(GCE). The deviations of u from unity in the hadron-resonance gas (HG) come due to Bose and 



ermi statistics, resonance decays (see, e.g., Ref. 



1^. 
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and exactly enforced conservations laws 



121]. Previous study [l2| has shown a rather rapid convergence of the scaled variances 
to asymptotic values. Finite system size effects can therefore be ignored for central collisions 
of heavy nuclei. 

The effect of Bose and Fermi statistics can be seen in the primordial (before resonance decay) 
values of the scaled variances in the GCE (see, e.g., Ref. ll|]). For chemical freeze-out at low 
collision energies (small temperature T and large baryonic chemical potential /i^) most charged 
particles are protons. Fermi statistics then lead to a suppression of positively charged and all 
charged particle number fluctuations, = Uch = 0.98 ^1. At high collision energies (large 
T and small ^b) most charged particles are pions. Thus, Bose statistics dominates and leads 
to the enhancement of particle number fluctuations, uj± = uJch = 1-05 ^ 1.06. These numbers 
demonstrate that the effects of quantum statistics are small at the chemical freeze-out. 

Resonance decay in the GCE and CE lead to the enhancement of particle number fluctua- 
tions. At high collision energies (the largest SPS energy and above) the GCE calculations give 
the following final (after resonance decay) values, uj± = 1.1 and Uch = 1.6 [11]. 

There is a qualitative difference in the properties of the mean multiplicity and the scaled 



2 



variance of multiplicity distribution in statistical models. In the case of the mean multiplic- 
ity results obtained with the GCE, canonical ensemble (CE), and micro-canonical ensemble 
(MCE) approach each other in the large volume limit (this reflects the thermodynamical equiv- 
alence of statistical ensembles). It was recently found [lO, 11, 12] that corresponding results 



for the scaled variance are different in different ensembles, and, thus, the scaled variance is 
sensitive to conservation laws obeyed by a statistical system. The differences are preserved in 
the thermodynamic limit. The global conservation laws imposed on each microscopic state of 
the statistical system lead to the suppression of particle number fluctuations. At large T and 
small fiB (high collision energies) the final state scaled variances behave in the CE as, uj±. = 0.8 
and uJch = 1-6 [ul, and in the MCE as, w± = 0.3 and oOch = 0.6 0. 

It has been found that the fluctuations in the number of nucleon participants give a large 
contribution to hadron multiplicity fluctuations in peripheral A+A collisions. Thus ideally a 
comparison between data and predictions of statistical models should be performed for results 
which correspond to collisions with a fixed number of nucleon participants. This can approx- 



imately be achieved by selecting only the most central A+A collisions (see Ref. 13|, [ij] for 
details). The conditions for the centrality selection in the study of fluctuations should thus 
be much more stringent than in the multiplicity measurements. To minimize the contribu- 



tions from the participant number fluctuations, one should restrict the analysis to the 



central A+A collisions selected by the number of projectile spectators (see Refs. [12 



% most 



15|). In 



this sample both the fluctuations of the number of participating nucleons and the resulting 
enhancement of multiplicity fluctuations should not exceed a few %. 

Usually the ideal HG gas is used for the particle multiplicity calculations. The extension of 
the ideal gas picture based on the van der Waals (VDW) excluded volume procedure have been 



suggested in Refs. [171, llSj, ll9|] to phenomenologically take into account repulsive interactions 



between hadrons. This leads to the suppression of hadron number densities and energy density. 
In the present paper we study the effect of the VDW excluded volume on the multiplicity 
fluctuations. The particle number ratios are independent of the proper volume parameter, if 
it is the same for all hadron species. Thus, a rescaling of the total volume V leads to exactly 
the same hadron yields as those in the ideal HG. In the present paper it is demonstrated that 
the multiplicity fluctuations are suppressed in the VDW gas. This suppression is qualitatively 
different from that of particle yields. In contrast to average multiplicities, the suppression of 
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multiplicity fluctuations can not be removed by rescaling of the total volume of the system. 

In the present study we restrict ourself to the GCE calculations within the Boltzmann 
statistics approximation. The paper is organized as follows. In Section II the ideal HG model 
is presented. In Section III the one-component VDW gas is used to calculate the average 
multiplicity and the scaled variance. These results are extended to the multi-component VDW 
gas in Section IV. Both primordial and final state (after resonance decay) hadron yields and 
their fluctuations are considered. In Section V the scaled variance of negative, positive and 
all charged hadrons is calculated along the chemical freeze-out line. This gives the VDW 
HG predictions for fluctuations of hadron multiplicities in central A+A collisions at different 
collision energies from SIS to LHC. A summary, presented in Section VI, closes the paper. 



II. IDEAL HADRON-RESONANCE GAS 



The partition function of the ideal Boltzmann HG in the GCE has the form: 



Hd 



exp 



where: 



A, 



27r2 



k'^dk exp 



+ m2)i/2 



T 



J 27r2 ^ ^ V T 



(2) 



In Eqs. ( ITpil . \j = exp {jij/T) is the particle specific fugacity, gj the degeneracy of jth particle 
species, rrij the particle mass, and iij = bjfiB + Sjfis + Qjfi'Q the chemical potential due to the 
jth particle internal quantum numbers {bj,Sj,qj). The j-summation in Eq. ([1]) is taken over 
all hadron species. Global chemical potentials (/i^, fis, fig) and temperature T regulate the 
systems average charges (baryonic, strangeness, electrical) and energy densities. Finally, V is 
the volume of the system and K2 is the modified Hankel function. 



A. Primordial Fluctuations 



In the GCE the first two moments of the multiplicity distribution of jth particle are: 



4 



Thus, (f)j is the jth hadron number density, n'j^. The variance of the jth particle number 
distribution is: 

{{AN,f) ^ {{N, - (iV,)))2 = (iV|) - {N,f = (iV,) , (5) 
and the scaled variance equals unity, 

_ ((AiV,)^) _ 

For the average multiplicity and variance of, for example, positively charged primordial hadrons 
(i.e. before resonance decay) one obtains: 

(N^) = E (^^■) ' ((^^+)') = E E i^N.AN,) . (7) 

j, qj>0 j, <jj>0 k, qk>0 

The results for the negatively charged or all charged hadrons are similar to Eq. (j7]) with summa- 
tion over negative charges, qj < 0, qk < 0, or all non-zero charges, qj 7^ 0, 7^ 0, respectively. 
For j k, one finds. 



and from Eqs. ([5]) and ([8]) it then follows 



d \ fd 



dXj /V dXk 



{AN.ANk) = {N,Nk) - (N,) (iV,) = 6,^ {N,) . (9) 

From Eq. iQ one concludes that in the ideal HG the correlations between different particle 
species are absent in the GCE. Any particle number distribution in the ideal Boltzmann gas is 
Poissonian (for large average multiplicity this is equivalent to a Gaussian) and its scaled variance 
therefore equals unity. For example, the scaled variance for positively charged hadrons is. 

This is also valid for negatively charged and all charged hadrons. 



B. Final State Fluctuations 

Final state yields and (co) variance have simple forms in the GCE [l^: 

{N,) = {n;) + J2 (Nr) ' (11) 

R 

{AN.ANk) = {{AN;Y)6,k + E ((AAT^)' + E {Nr ){An,Ank)R , (12) 



R R 
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where N* is primordial j'th hadron multiplicity, and: 

{n,)R = J2^r , (13) 

r 

{AnjAnk)ji = {njnk)R - {nj)R{nk)R = ^ fef rf-^^ nl j^ - {nj)ji{nk)R , (14) 

r 

is the branching ratio of the rth decay channel of Rih resonance with normalization condition 
Ylir^r^ ~ 1' while rfj j^ and n], R are the respective multiplicities of particles of species j and k 
in rth decay channel of i?th resonance. Generally, one finds enhancement of fluctuations due 



to resonance decay which leads to ujj > 1 (valid in GCE and CE, see, e.g., Ref. [ll|). The first 
term in the right hand side of Eq. (fT2|) corresponds to fluctuation of i?th resonance multiplicity 
and the second term appears due to the probabilistic character of resonance decay. For the 
Boltzmann gas discussed in this paper, it follows, {{ANrY) = {Nr), {{ANj)"^) = {N*), and 
Eq. ( fT2|) simplifies to: 



{AN.ANk) = {N;)6,k + 2^ {NR){n,nk)R. (15) 

R 



From Eqs. (JTTl) and ( 1151) it immediately follows, 

^ (N;) + Er {Nr ){n])R 

(N;) + Er {NR){n,)R ■ ^''^ 

Thus, Ljj > 1 as a consequence of the inequality (^j,i?)^ > n^j^^^ , i-e. the scaled variance of jth 
hadron exceeds 1 due to the presence of resonances decaying into more than one jth hadron. 
Using Eqs. (fTTl) and ([T5|) for average yields and correlators one finds from (JT]), 

^ E„,,>o {{N;) + J:,,,,>oJ:r {NR){n,n,)R] 

E,,,,>o [(N;) + Er {NR){n,)R] ' ^ > 

and similar expressions for final state cj. and Uch with the summation over < 0, qk < and 
7^ 0, qk ^ 0, respectively. 



III. ONE-COMPONENT VDW GAS 

The VDW excluded volume procedure gives the GCE partition function of one component 
Boltzmann gas [17\: 

Z(V^,T,/i,) = ^ — [(V - v,N) 0, f eiV-v,N) , (18) 

7V=0 



6 



where vj is the proper volume of the jth particle. For the 'hard sphere' particles with radius rj 
the volume parameter vj equals the 'hard-core particle volume', 47rr|/3, multiplied by a factor 
of 4. A Laplace transform of Eq. f|T8|) reads: 



Z{s,T,fij) = / dVexp{-sV) Z{V,T,fij) 
'o 



N=0 



dVexp{-sV) {V - VjN) 



N 



ViN 



exp[—VjS 



The system pressure is defined by the pole- singularity, s*, of the function, Z (fTOl) . 

\nZ 



p = lim T = T s* 

V^oo V 



and, thus, satisfies the following equation, 

p = exp 



(19) 



(20) 



(21) 



For Vj = 0, Eq. (l2T]) is evidently reduced to the ideal gas result, pid = Tcpj. The average particle 
number in the VDW gas equals to: 

To find the particle number fiuctuations one needs to calculate. 



(22) 



^'^ - -z [^'k 



dp 



V 

T ^ d\] 



(23) 



The variance is therefore: 



dp 



+ A 



d'^p 



(24) 



Thus, one can write the scaled variance as the following. 



From Eq. fl?I]) one finds: 



dp 



dp 
dX] 



-1 



(25) 



A 



dX, 



T Xi 



dp_ 



where the notation, 



Xj = exp ( — 



Vj p 



'3 1 



(26) 



(27) 
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is introduced to simplify the formulas below. From Eq. (l26i) it follows, 

dp T Xj 



A. ^ = , : ^ , (28) 
oXj 1 + VjXj 



and one finds: 



From Eqs. ( l2T]|27|29p one finds the VDW equation of state, familiar from textbooks, 

P[V - v,{N,)] = (N,) T . (30) 
At small particle density, vjXj <^ 1, one finds from Eq. (l29ll . 

{N,) = V 0, (1 - 2v,<l>,) = (iVf ) (1 - 2v,r^f) , (31) 

where n^J^ = is the ideal gas particle number density. In the opposite limiting case of high 
particle density, VjXj ^ 1, one obtains, 

iV, = - , (32) 

thus, the value of Ijvj is the upper limit for the particle number density in the VDW gas ( |T8l) . 
From Eq. ([26]) and ([28]) it follows. 



A" 



2 d^p _ Tv,x] 



' d\] (1 + VjX^Y 



2 _ ^i^i 



1 + VjXj 



(33) 



From Eqs. fl25|28ll33p one then obtains, 

^3 = (1 + ^jV^ ■ (34) 

From Eq. ( 134]) it follows that the scaled variance in the VDW gas is always smaller than the 
ideal gas value ([6]) ujj = 1, i.e. the non-zero proper volume Vj suppresses the particle number 
fluctuations. At small particle number densities, VjXj <^ 1, the particle number fluctuations 
in the VDW gas are approximately equal to those in an ideal gas IQ. The Eq. (IMl) gives, 
ujj = 1 — '^'^j'Pj ■ Oil the other hand, at VjXj ^ 1 the scaled variance of the VDW gas goes 
to zero as Uj = {vjXj)""^. In this limit the particle number density is close to its lower limit, 
rij = Thus, there is no 'free space' for particle number fluctuations. 
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IV. MULTI-COMPONENT VDW GAS 



The partition function of the multi-component VDW gas equals to 



17|,118|,119|: 



CO _ 

j=l Nj=0 ^ 



0{y - I^^^.AT,) . 



(35) 



Similar to Eq. (l2Ti) one finds the equation for the pressure function in the multi-component 
VDW gas, 



P = Y. T 



X 



3 ' 



where Xj is given by Eq. ([2 



(36) 



A. Primordial Fluctuations 



The average multiplicity of jth particle is: 



_d_ 
dXi 



InZ 



V dp 
T 'dX, 



Similar to Eq. ( 12^ one obtains: 



(AiV,AiV,) = ^ 



dp d'^p 



dX, 



dXjdXk 



Calculating the derivative from Eq. ([36 

dp 



X, 



T Xj 



dX, 



and using Eq. ( 1371) one finds: 



1 + Y.i Vi Xi 

Vx, 



(37) 



(38) 



(39) 



(40) 



1 + Yji ^i^i' 

which extends Eq. f l29p to the multi-component VDW gas. Using Eq. fl39l) for the first derivatives 
of the pressure one finds: 



XA 



d'^p 



T XjXj^ 



Vj + Vk 



dX.dXu {l + T.iXiv{) 

The Eq. ( HTj) for j = k leads to the following result for Uj, 



(41) 



Uj = 1 + Xj 



d'^p 
9Af 



X 



dp \ ^ 



dX, 



1 - 



1 + 



1 + ^iViXi 



■ (42) 
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The Eq. (H2l) is reduced to Eq. ( l34l) for the one-component VDW gas. For an 'almost' ideal gas 
when all proper volumes come to zero, f j — 0, one finds the ideal gas results, {Nj) = V(j)j and 
LVj = 1, from Eqs. fHOl) and fj42l) . respectively. The above equations are simplified if the proper 
volumes of different hadron species are equal to each other, Vi = v. The scaled variance (142!) 
becomes then equal to: 

_ Xj {2 + V 

^ " (1 + V ■ ^''^ 

The Eq. ( H3|) demonstrates a suppression of the particle number fluctuations for all hadron 
species {uj < 1). At small total particle number density, v'^-Xi <^ 1, it behaves as, = 
1 — 2v(f)j, similar to the one-component VDW gas. In the opposite limiting case, vY^-Xi ^ 1, 
one flnds, Uj = 1 — (f^j/J^i'Pi = 1 ~ ''^f / ''^'tot- ^ dense VDW gas the particle number 
fluctuations are suppressed, and the suppression for jth species is proportional to the ratio of 



the ideal gas jth particle number density, n^f", to the total particle number density, nl'^^ = nf. 
The largest suppression takes place for the total multiplicity, Utot — ^ at t> Xj oo. This is 
in an agreement with the behavior of the scaled variance in one-component VDW gas at high 
particle number density. 

For the correlators (|38l) in the VDW gas one obtains. 



'jk 



Xk (2 + vYj^Xi 



(44) 



(1 + ^E.^.)' 

In contrast to the ideal HG result ([9]), the excluded volume effects lead to the (anti) correlation 
between different particle species, j ^ k, in the VDW gas seen from Eq. fj44|) . The physical 
origin of these anticorrelations is rather clear. The 'large' number of jth particles, ANj > 0, 
reduces the available free space for kth particles. This makes preferable to have a 'small' number 
of kth. particles, ANk < 0, thus, leading to {ANjANk) < 0. 

Using Eq. (jUj) one flnds the primordial fluctuations of the positively charged hadrons, 

E,, ,,>o E., ,,>o (AAT, AiV,) 



E 



i, qi>0 



1 



E E 

i,(?j>Ofc,gfc>0 



Xj Xk{2 + V Xj 
(1 + vJ2^Xif 



X 



E 



(45) 



and similar expressions for uj_ and Uch with the summation over qj < 0, < and qj ^ 
0, Qk 7^ 0, respectively. 
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B. Final State Fluctuations 



To take into account resonance decay we follow the procedure from Refs. ll|, ll2|. The corre- 
lators for the final hadrons are expressed as, 

(ANjANk) = {AN* an;) + J](iV^) K-n,)^ + ^^(AiV; AiV^) K)^ 

R R 

+ J2{AN; ANn) {n,)R + (^^R ^^R') (^:i)R ■ (^6) 

R R^R' 

where all particle-particle, (ANjAN^), particle-resonance, {AN*ANj^), and resonance- 
resonance, (ANjiANri) , correlators in Eq. fH6|) are calculated according to Eq. (jUj). Note 
the essential difference between the correlators in the ideal gas (ITSl) and those in the VDW gas 
(H6l) . The new terms in Eq. ( l46i) come from the anticorrelations of different stable hadron and 
resonance species in the VDW gas according to Eq. (jH 



V. SCALED VARIANCE ALONG CHEMICAL FREEZE-OUT LINE 



In this section we present the results of the VDW HG for the scaled variances along the 



chemical freeze-out line in central A+A collisions for the whole energy range from SIS to 



The procedure to define the chemical freeze-out line is essentially the same as in Refs. 



11 



HC. 



12| 



The values of T and /i^ at the chemical freeze-out at different collision energies are presented 



in Table I. They are almost identical to those values in Fig. 1 and Table I of Ref. 



ll|. The only 



tiny difference comes because of the Boltzmann statistics approximation used in the present 



20| 



paper. Note that the conditions for average energy per particle, {E) / {N) = 1 GeV 
zero value of the net total straiigeness, 5 = 0, and the charge to baryon ratio, Q/B = 0.4 



remain the same as in Refs. 



11 



12|. This procedure is possible because all particle ratios 



and energy to particle ratio in the VDW HG gas (with the same hard-core volume Vj 



for all hadron species) remain unchanged in a comparison to the ideal HG, fij/uk = n^f/n^ 
and e/ritot = ^^'^/^m- The dependence of on the collision energy is parameterized as ij. 
H-B (a/^7V7v) = 1-308 GeV (1 + 0.273 ys/vw)^^ , where the cm. nucleon-nucleon collision energy. 
^/sNN, is taken in GeV units. The strangeness saturation factor, 75, is parameterized as [2| 
75 = 1 — 0.396 exp (— 1.23 T/fis)- Both these relations are the same as in Refs. tU 



El 



The excluded volume corrections, i.e. the factors exp {—vp/T) and (1 + vY2i 

cal- 



culated using the THERMITS package 



2l|. Numerical optimization functions allow to solve 
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\/snn 


T 




is 


R = 




rii/ntot 


[GeV] 


[ MeV ] 


[ MeV ] 




r =0.3fm 


r =0.5fm 






nch/ntot 


2.32 


64.3 


800.8 


0.64 


0.944 


0.788 


n /1 1 Q 
U.4io 
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n A fi7 
U.4D ( 


4.86 


116.5 
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0.870 
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U.ooo 


U.ioD 


U.D ( o 


6.27 


128.5 


482.4 


0.72 


0.844 


0.552 


U.o ( i 


U.ZU ( 


n 1^77 

U.D ( ( 


7.62 


136.1 


424.6 


0.74 


0.825 


0.519 


U.ooo 


U.zio 


u.o (0 


8.77 


140.6 


385.4 


0.75 


0.812 


0.498 


U.o4y 


U.zzo 


u.o / 4 


12.3 


149.0 


300.1 


0.79 


0.786 


0.459 


U.ooi 


U.Zo ( 
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17.3 


154.4 


228.6 


0.83 


0.766 


0.432 


0.316 


0.245 


0.561 


62.4 


160.6 


72.7 


0.98 


0.738 


0.397 


0.285 


0.263 


0.548 


130 


161.0 


35.8 


1.0 


0.735 


0.393 


0.278 


0.268 


0.546 


200 


161.1 


23.5 


1.0 


0.735 


0.393 


0.277 


0.269 


0.546 


5500 


161.2 


0.9 


1.0 


0.735 


0.393 


0.273 


0.273 


0.546 



TABLE I: The chemical freeze-out parameters T, fiB, and 75 in central A+A collisions are presented 
at different cm. energies y^SAw- The calculations are done for the HG with Boltzmann statistics 
in the GCE. The excluded volume parameter, v = 167rr'^/3, is taken to be the same for all hadron 
species. The VDW suppression factor, R = exp{—vp/T)[l + vY^-Xi], for the VDW energy density 
and particle number densities, and the ratios rii/ntot = ''T'l'^/iT'tot (^^^ ^^S- 2, right) are presented. 

the transcendental equation fl2T|) for the pressure and, thus, find the hadron yields and fiuctua- 
tions of the VDW HG. The thermodynamical parameters T, fis, and 75, the VDW suppression 
factor R = exp(— t>p/T)[l + v Xj], and the particle number ratios rii/ritot = /nl'^^ (see also 
Fig. 2, right) along the chemical freeze-out line are presented in Table I. Once a suitable set 
of thermodynamical parameters is determined for each collision energy, the scaled variance of 
negatively, positively, and all charged particles can be calculated using Eq. ( H5l) for 0;+ and 
similar equations for uj_, Uch, with the correlators taken from Eq. fj46l) . The resulting scaled 
variances are presented in Tables II-IV and Figs. 1 and 2 (left) as a function of ^/sNN. The 
values of ^/snn quoted in Tables I-IV and marked in Figs. 1-2 correspond to the beam energies 
at SIS {2 A GeV), AGS (11. 6 A GeV), SPS {20 A, 30 A, 40 A, 80 A, and 158 A GeV), colliding en- 
ergies at RHIC {y^sj^ = 62 A GeV, 130 GeV, and 200 GeV), and LHC {^/sj^ = 5500 GeV). 
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FIG. 1: The scaled variances for positively {left) and negatively {right) charged particles along the 
chemical freeze-out line for central A+A collisions (see also Tables II and III). The calculations are 
done for the HG with Boltzmann statistics in the GCE. Symbols at the lines correspond to the specific 
collision energies pointed out in Table I. The dash-dotted lines correspond to the primordial and solid 
lines to final state values. The pair of upper lines show the ideal HG results, the pair of middle 
lines the VDW HG with hard-core radius r = 0.5 fm. The lowest pair of lines corresponds to largest 
possible VDW suppression for the primordial (dotted) and final (solid) scaled variance, respectively. 
This largest VDW suppression happens at ntot = 1/v. 

Both the primordial and final state scaled variances are presented. To make a correspondence 
with real measurements, both strong and electromagnetic decays should be taken into account, 
while weak decays should be omitted. 

Some features of the results should be mentioned. The values uj± = Uch = 1 are shown by 
the dashed-dotted lines in Figs. 1-2. They correspond to the primordial fluctuations in the 
ideal HG (r = 0). The bump in for final state particles seen in Fig. [T]at the small collision 
energies is due to a correlated production of proton and vr"*" meson from decays. This 
single resonance contribution dominates in at small collision energies (small temperatures), 
but becomes relatively unimportant at the high collision energies. 

The number of negative particles is relatively small, (A^_) <^ (^+)) at low collision energies 
(see Fig. 2, right). For the hard-core radius r = 0.5 fm, this region corresponds 'small' particle 
number density, vritot <^ 1. It then follows, for the primordial values, uj± = 1 — 2vn± and 
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1.183 


1.106 


0.970 


0.681 


7.62 


1 


0.939 


U.oo / 


n RAO 


1.178 


1.091 


0.948 


0.671 


8.77 


1 


0.936 


U.OOO 


U.oOi 


1.173 


1.079 


0.932 


0.663 


12.3 


1 


0.931 


U.ooZ 


u.Doy 


1.158 


1.051 


0.898 


0.645 


17.3 


1 


0.928 


0.832 


0.684 


1.145 


1.029 


0.873 


0.632 


62.4 


1 


0.928 


0.840 


0.715 


1.115 


0.989 


0.835 


0.617 


130 


1 


0.929 


0.843 


0.722 


1.109 


0.983 


0.830 


0.617 


200 


1 


0.929 


0.844 


0.723 


1.107 


0.981 


0.830 


0.617 


5500 


1 


0.930 


0.845 


0.727 


1.103 


0.978 


0.828 


0.617 



TABLE II: The scaled variances of the primordial and final positively charged hadrons in central A+A 
collisions are presented at diff'erent cm. energies ^/SNN (see also Fig. 1, left). The calculations are done 
for the HG with Boltzmann statistics in the GCE. The excluded volume parameter, v = 167rr^/3, is 
taken to be the same for all hadron species (r = corresponds to the ideal HG results). The condition 
ntot — 1/u corresponds to largest possible VDW total particle density, and this gives an upper limit 
of the VDW suppression effect for the scaled variances. 

LOch = 1 — ^vrich- Thus, the VDW suppression effects are seen at low collision energy for 
positively charged and all charged particles (Fig. 1 and 2), and are absent for negatively charged 
particles (Fig. 1, right). 

At high collision energies (large T) one finds n_ = n+ as seen from Fig. 2, right. For the 
hard-core radius r = 0.5 fm. Fig. 1 shows the asymptotic values of 0.85 for the primordial and 
0.83 for the final scaled variances u;± in the VDW HG, instead of 1 and 1.1, respectively, in 
the ideal HG. Thus, one observes about 15% and 25% VDW excluded volume suppression for, 
respectively, the primordial and final values of cj-(-. The upper limit of the VDW suppression 
effects for the scaled variances is obtained at rifot — ^/v and is presented in Figs. 1-2. The 
approximate relations both at small density, a; j = 1 — 2vnj, and high density, a;j = 1 — rij/ntot, 
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Primordial 






Final 




[GeV] 


r = 


r =0.3fm 


r =0.5fm 


ntot = V ^ 


r = 


r =0.3fm 


r =0.5fm 


ntot = V ^ 


2.32 


1 


0.997 


U.yoy 


U.y4D 


1.001 


0.998 


0.988 


0.941 


4.86 


1 


0.976 


u.you 


naif; 

U.oiO 


1.016 


0.984 


0.920 


0.763 


6.27 


1 


0.968 


019 

u.yiz 


n 70^ 
U. Ivo 


1.026 


0.980 


0.898 


0.723 


7.62 


1 


0.963 


U.yUi 


n VQO 
U. loz 


1.035 


0.977 


0.883 


0.700 


8.77 


1 


0.957 


u.<5y4 


U. M 


1.041 


0.976 


0.873 


0.686 


12.3 


1 


0.951 


U.ooU 


U. (Do 


1.056 


0.973 


0.855 


0.660 


17.3 


1 


0.944 


0.870 


0.755 


1.068 


0.972 


0.843 


0.644 


62.4 


1 


0.933 


0.852 


0.737 


1.092 


0.973 


0.827 


0.621 


130 


1 


0.931 


0.849 


0.732 


1.097 


0.975 


0.827 


0.619 


200 


1 


0.931 


0.848 


0.731 


1.099 


0.976 


0.827 


0.618 


5500 


1 


0.930 


0.846 


0.727 


1.103 


0.978 


0.828 


0.617 



TABLE III: The same as in Table II, but for negatively charged hadrons (see also Fig. 2, right). 




FIG. 2: Left. The same as in Fig. 1, but for Wc/j (see also Table IV). Right. The ratios Ui/ntot = n^f / n^tot 
(see also Table I) for i = + (dashed line), i = — (dotted line), and i = ch (solid line). 

demonstrate that VDW suppression of jth fluctuations is proportional to jth density in the HG. 
This explains why the VDW suppression for the primordial ujch seen in Fig. 2 is approximately 
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Primordial 






Final 




[GeV] 


r = 


r =0.3fm 


r =0.5fm 


ntot = V ^ 


r = 


r =0.3fm 


r =0.5fm 


ntot = V ^ 


2.32 


1 


0.974 


u.yu4 


U.Ooo 


1.061 


1.034 


0.960 


0.574 


4.86 


1 


0.927 


U. ( oo 


n /I97 


1.331 


1.236 


1.049 


0.585 


6.27 


1 


0.912 


U. ( 00 


U.4ZO 


1.397 


1.274 


1.056 


0.592 


7.62 


1 


0.901 


U. / oo 


n /1 0/1 

U.4z4 


1.440 


1.296 


1.058 


0.598 


8.77 


1 


0.895 


U. / zo 


U.4ZD 


1.468 


1.309 


1.059 


0.602 


12.3 


1 


0.882 


n 71 9 


U.4oZ 


1.521 


1.331 


1.059 


0.612 


17.3 


1 


0.873 


0.702 


0.439 


1.557 


1.344 


1.059 


0.619 


62.4 


1 


0.861 


0.692 


0.452 


1.601 


1.355 


1.056 


0.632 


130 


1 


0.860 


0.691 


0.454 


1.604 


1.356 


1.056 


0.634 


200 


1 


0.860 


0.691 


0.454 


1.605 


1.356 


1.056 


0.634 


5500 


1 


0.860 


0.691 


0.454 


1.605 


1.356 


1.056 


0.634 



TABLE IV: The same as in Table II, but for all charged hadrons (see also Fig. 2, left). 



the same as for a;+ at small collision energy and 2 times larger than that for uj± at high collision 
energies. 

A comparison of the primordial scaled variances with those for final hadrons demonstrates 
that the fluctuations increase in GCE ideal HG for each hadron species (see Eq. f ll6p ) for 
positive, as well as for negative, and all charged hadrons (see Figs. 1-2). In the VDW HG the 
behavior is more complicated. At high collision energies, as demonstrated in Fig. 1, they even 
decrease for uj± after resonance decays. This unexpected behavior follows from Eqs. (jH]) and 
( 146|1 . Calculating fH6l) one finds that only one term, '^ji{Nfi){ninj), coming from resonance 
decays gives a positive contribution to the (ANjANk) correlator. Other terms in the right hand 
side of Eq. ( H6l) appear due to particle- resonance and resonance-resonance (anti) correlations in 
the VDW HG. As seen from Eq. these terms are negative and, thus, suppress the scaled 
variances for the final state particles. The resulting effect from resonance decays is defined 
by the competition of two effects: decays into more than one jth particle against particle- 
resonance and resonance-resonance anticorrelations. Multi-particle decays lead to increase the 
fluctuations, while the anticorrelations because of the excluded volume lead to decrease them. 
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For all charged particles the multi-particle decays win and ooch always increases due to resonance 
decays as seen from Fig. 2. This same is true for a;+ at small collision energies. The opposite 
situation takes place for ui± at high collision energies. Anticorrelations because of the excluded 
volume overcome the multi-particle decay contributions, and u!± decreases as seen from Fig. 1. 

VI. SUMMARY 

The hadron multiplicity fluctuations in relativistic nucleus-nucleus collisions have been con- 
sidered in the statistical hadron-resonance gas model. We study the effect of the van der Waals 
excluded volume on the hadron distribution scaled variances. In the present paper we restrict 
ourself to the grand canonical calculations within the Boltzmann statistics approximation. 

If the proper volume parameter is the same for different hadron species, the particle number 
ratios remain unchanged and equal to those in the ideal hadron-resonance gas. Therefore, with 
rescaling of the total volume V one may compensate the excluded volume suppression of the 
hadron densities leaving the hadron yields the same as in the ideal hadron-resonance gas. In 
the present paper it has been demonstrated that the multiplicity fluctuations arc suppressed 
in the van der Waals gas. This suppression is qualitatively different from that of the particle 
yields. In contrast to the average multiplicities, the suppression of multiplicity fluctuations can 
not be removed by rescaling of the total volume of the system. 

In this work we have considered two 'reasonable' values of hard-sphere radii, r = 0.3 fm and 
r = 0.5/m, as well as the limiting behavior of the VDW suppression. Estimates of possible 
excluded volume suppression effects of particle number fluctuations from existing data in A-l-A 
collisions will be the subject of future studies. 
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